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Using Monte Carlo simulations, we study a Higgs transition in several three-dimensional lattice 
realizations of the noncompact CP 1 model (NCCP 1 ), a gauge theory with two complex matter fields 
with SU(2) invariance. By comparing with a one-component theory, which is well understood and 
has continuous transition in the inverted XY universality class, we argue that the two-component 
case also has continuous Higgs transition with a larger correlation length exponent (i.e., it is "more 
continuous"). The transition can become first order in the vicinity of a new "Molecular" phase, 
which occurs in one of our models, but is continuous in a wide range of parameters away from this 
phase. The situation is significantly clarified by studying a model where the Molecular phase is 
entirely absent, and a wide regime with a continuous transition can be readily established. The 
two-component theory is also an effective description of the hedgehog-suppressed 0(3) universality, 
and results are relevant for the recently discussed "deconfined quantum criticality" scenario for the 
continuous Valence Bond Solid to Neel quantum phase transition. 



I. INTRODUCTION 

In this paper, we examine the nature of the Higgs tran- 
sition in a two- component theory with a dynamical gauge 
field in three (classical) dimensions (3D)ii A schematic 
continuum action is 



(V-m)*| 2 + m|*| 2 + u|*| 4 + K (Vxa) 2 . (1) 



^ has two complex components and a is the gauge field 
that can take any real values (i.e., it is non-compact). 
This theory is also called the noncompact CP 1 model 
(NCCP 1 )^ When the charged field * is gapped, the 
gauge field is massless and there is a free photon, but 
when \& condenses, a becomes massive - this is the Higgs 
transition. In Ref. 0, we argued that this action also 
describes an 0(3) classical spin system with prohibited 
hedgehog topological defects on long scales. In particu- 
lar, the ordering transition in the hedgehog-suppressed 
model and the above Higgs transition are in the same 
universality. The nature of the transition is of interest 
in the context of the recently discussed quantum phase 
transitions that lie outside the Landau-Ginzburg- Wilson 
frameworki&ii^i&i Thus, it has been proposed^^ that 
the Neel to Valence Bond Solid transition in a spin-1/2 
system with SU(2) symmetry can undergo a continuous 
quantum phase transition described by the same univer- 
sality Eq. H|). 

Here we present a detailed numerical analysis of sev- 
eral lattice realizations of the theory ((T|) (we will refer to 
these as lattice superconductor models^) . We extend the 
study of Ref. of the SU(2)-symmetric two-component 
model, with the goal to check more thoroughly the pro- 
posed second-order nature of the transition. Here, we 
argue that our models have continuous Higgs transition 



over wide parameter ranges, see Figs. QJ>,c. Note that 
in one of our models, Fig. [TJd, we also find a parameter 
range, in the vicinity of a third "Molecular" phase [pro- 
duced when gauge neutral pairs of the Vf^ particles of 
Eq. ((T|) condense] , where the transition becomes first or- 
der. This is always a possibility, and in the specific model 
happens generically near this new phase. We adopt two 
fresh strategies to establish the nature of the transition 
in these systems. First, in the regime where we believe 
the transition is second-order, we compare with a one- 
component model for which it is well established that 
the transition is continuous^ and show that the tran- 
sition in our two-component models is "more continu- 
ous", in the sense that the thermodynamic singularities 
are more weak. Second, we introduce a lattice realiza- 
tion of the model in Eq. |T]), where the Molecular phase 
is entirely absent. The Higgs transition is then shown 
to be continuous over the entire phase boundary, which 
greatly simplifies the interpretation. Our best estimates 
for the exponents at this transition are v — 0.7 — 0.75 
and rj = 0.2 — 0.4. Recent Quantum Monte Carlo simula- 
tions of square lattice S = 1/2 Hciscnbcrg antifcrromag- 
nets with four spin exchange terms and SU(2) symmetry 
have also found a continuous relativistic transition be- 
tween a Neel and Valence Bond Solid phased which 
is believed to be in the same universality class as the 
transition studied here. The exponents obtained in those 
studies arc r/ = 0.26(3) and v = 0.78(3) in Ref. i and 
j] = 0.35(3), v = 0.68(4) in Ref. i, consistent with the 
ones reported here. 

Other studies,— ' 10 i n particularly in two-component su- 
perconductor models where the SU(2) symmetry is bro- 
ken down to U(l) (easy plane NCCP 1 ), have suggested 
that such Higgs transitions may always be first order. 
While those authors so far have found first order tran- 
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sitions in these U(l)xU(l) models, we do not know any 
fundamental reasons why this needs to hold generally. 
This is particularly true in light of the evidence presented 
here for a continuous transition in the S U (2) symmetric 
model. Moreover, the studies in Refs. Hn] were per- 
formed on models where the Molecular phase is present, 
which we suspect may be driving the first order nature of 
the transition. In future work we will apply the ideas de- 
scribed here to obtain and study models without a Molec- 
ular phase in this symmetry class as well, which should 
help clarify the issue. 

The study of the criticality in the theory JT]) dates 
back to Halperin, Lubensky, and Mai At the fixed 
point where the matter and gauge fields are decoupled 
(e 2 ~ 1/k — > 0), the gauge interaction is a relevant per- 
turbation, and no new stable fixed points were found in 
d — 4 — e treatment when the number of components N is 
less than 365. This was interpreted as an evidence for a 
fluctuation-induced first-order transition. For larger TV, 
and also in a l/N treatment directly in d = 3, the transi- 
tion was found to be continuous. Later work by Dasgupta 
and Halperin^ showed that the transition in the N = 1 
case can also be continuous. The one-component super- 
conductor model is dual to the usual 3D XY model that 
has only short-range interactions. Concisely stated, one 
can fruitfully think about the one-component system in 
terms of the Abrikosov-Nielsen-Olesen (ANO) vortices 
which carry quantized gauge flux and have only short- 
range interactions. The duality claim then follows if we 
recall that the 3D XY model can be also viewed as a sys- 
tem of short-range interacting current loops. The Higgs 
phase where the ANO vortices are gapped corresponds to 
the disordered phase of the XY model with small loops, 
while the disordered phase of the charged superconduc- 
tor where the ANO vortices condense corresponds to the 
ordered phase of the XY model with proliferated loops. 
Thus, the ANO vortices play a crucial role, which is dif- 
ficult to capture in the 4 — e treatment. We expect that 
the same physics happens in the two-component case, 
but with the ANO vortices that now have non-trivial in- 
ternal structure. While our analytical picture is not as 
complete as in the one-component case, since we are un- 
able to perform a duality transformation in the present 
case, one of our models has a limit that strongly suggests 
such an interpretation, and we confirm that the transi- 
tion is continuous by direct Monte Carlo studies on the 
lattice. 

One of our main strategy exploring such new two- 
component matter - gauge systems in Monte Carlo is to 
compare closely with the well-understood one-component 
case. We simulate both systems under similar conditions 
and compare directly obscrvablcs in the gauge sector that 
detect the Meissner-Anderson-Higgs physics of charged 
condensates. We also study separately new features in 
the two-component system, such as the presence and pos- 
sible condensation of composite gauge-neutral fields. 

The paper is organized as follows. In Sec.HU we define 
the models and describe their overall phase diagrams. In 



Sec. IIII1 we focus on the Higgs transitions and present 
several two-component cases where we argue the tran- 
sition is continuous. We also discuss the observed sys- 
tematics along the full Photon to Higgs phase boundary 
and show that we can improve it by suppressing possible 
other phases. In Sec. HV[ we summarize the results and 
conclude with an outlook. 



II. MODELS AND THEIR PHASE DIAGRAMS 

We study the following three models. The first model 
realizes a one-component lattice superconductor: 

Snccp° = -J"^cos(0j - l+jl + Oi M ) + S a , (2) 
i . fj 

S a = ^(Axa) 2 . (3) 

□ 

A U(l) matter field e resides on the sites of a 3D cu- 
bic lattice, while a non-compact gauge field a M resides on 
the links; S a is the usual lattice Maxwell action for the 
gauge field. We will use 'NCCP ' as a short-hand refer- 
ence to this model, where 'CP ' refers to one matter field 
component while 'NC emphasizes that the gauge field is 
non-compact.^ This model and the ones below have two 
parameters J and K,. The NCCP model is our reference 
system for matter-gauge simulations. 

The second model realizes a two-component lattice su- 
perconductor: 

SWcCP 1 , model I = { Z \ Z i+fi elai * + c - c ) + S " ■ 

(4) 

Here the matter field is a CP 1 field z = (z-f,zi), \z\ = 
1 . This model was introduced in Ref . [H and studied at 
JC = 0.6 using small system sizes; the present work is an 
extension of that study. 

The third model is a different realization of the two- 
component case: 

•SWcCP 1 , model II = ^NCCP 1 , model I + ^ m ^0 (J\ z i z i+fi I ) ■ 

(5) 

Here Io(x) is the modified Bessel function. Compared 
to the model I, there is an additional antiferromagnetic 
interaction between nearest-neighbor spins. This interac- 
tion eliminates the unwanted 'Molecular' phase from our 
phase diagram. Our motivation for this model will be- 
come clear soon, while here we only note that the specific 
form allows some analytical understanding of the phase 
diagram and that any such short-range interaction that 
respects the symmetries of the problem does not change 
the universality of the transition. 

Figure [1] shows the phase diagrams of the three models 
in the K, — J parameter space. In each case, for small 
J , the matter fields are gapped and the gauge field is 
masslcss described by the Maxwell term - we mark this 
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as the Photon phase. On the other hand, when J and 
JC are sufficiently large, the matter fields condense and 
gap out the photon - this is the Higgs phase. In the one- 
component case, there remain no gapless excitations. In 
the two-component system, we can construct a gauge- 
invariant 0(3) vector field 



n = z^az . 



(6) 



This obtains an expectation value in the Higgs phase, and 
there appear Goldstone modes associated with the global 
symmetry breaking. The focus in this paper is on the 
Photon to Higgs transition driven by the condensation of 
the charged matter fields. We will also see how to view 
this transition coming from the Higgs side as driven by 
proliferation of the ANO vortices with internal structure. 

Proceeding with more details on the phase diagrams, 
consider first the NCCP model. At K = oo, there is no 
gauge field left and the transition along this axis is 3D XY 
ordering of the e l9 field (in Fig. QJi, this transition point 
is indicated with a star symbol). When JC becomes finite 
and the matter-gauge coupling is switched on, the univer- 
sality of the transition changes. Dasgupta and Halperin* 
argued that sufficiently far from the JC = oo limit, the 
transition becomes inverted XY as summarized in Sec. U 
in terms of the ANO vortices. We are also mostly in- 
terested in moderate JC. Using classical Monte Carlo, 
we studied several points along the phase boundary and 
found continuous transitions in our model; these points 
are indicated with open symbols. (The limit of large but 
finite JC was not answered in Ref. |8|, but several scenarios 
were discussed. One possibility is that the whole line is 
described by the inverted XY universality.) 

The picture of the transition becomes particularly sim- 
ple in the limit J = oo. Here, the gauge field takes on 
2-7r-integer values and the model reduces to that of con- 
served discrete- valued fluxes = (V x a) M with steric 
interactions, 



^NCCP 



E 



exp 



B=27rxint,V B=Q,B t , 



■w > B 

□ 



(7) 

This is the familiar loop representation of the 3D XY 
model and has been studied in great detail in the litera- 
ture. Note that in our ensemble derived from the gauge 
model the total flux is constrained to be zero, but this 
is not important in the thermodynamic limit. There are 
only small loops for JC > JC C = 0.07606, but they prolif- 
erate for JC < JC C . The transition is 3D XY in terms of 
these short-range interacting loops, which are precisely 
the ANO vortices carrying quantized flux. These loops 
are gapped on the Higgs side and condense upon reduc- 
ing JC and J , which is why the transition is referred to 
as 'inverted XY' (as opposed to 'direct XY' transition at 
JC = oo where currents conjugate to e proliferate upon 
increasing J) . The universality of the transition remains 
unchanged also away from the J = oo limit. Indeed, one 
can see that finite J only leads to short-range modifica- 
tions of the ANO vortex interactions and therefore does 



not change the nature of the transition. To conclude the 
description of the NCCP phase diagram, we readily see 
that there are no transitions along the JC = or J = 
axes in Fig. QJi. 

Consider now the NCCP 1 model I, whose phase dia- 
gram is shown in Fig. [TJd. Similarly to the NCCP case, 
there is no transition along the J = axis. Along the 
JC = oo line, the gauge field is again absent and the transi- 
tion is in the 0(4) universality class, but the universality 
changes once JC becomes finite. The region where we find 
the transition is continuous in a new universality class is 
marked with open symbols in the figure (details of the 
numerical study of the transitions are given in Sec. IIII[) . 

Unlike the one-component case, this two-component 
model has an additional phase at small JC and large J7, 
which, while interesting on its own, does not bear on 
the universality class of the Higgs transition that we are 
interested in. For completeness and in order to better 
understand the numerics in this model, let us describe 
the new phase, which we mark as Molecular. Here, the 
charged field z is gapped but the charge-neutral field n 
obtains an expectation value, and the 0(3) order and 
the photon coexist. This molecular condensate appears 
because of the attraction between oppositely charged z^ 
and z fields. Explicitly, consider the JC — axis. The 
gauge field can be integrated out completely, with the 
resulting action for the z's 



^NCCP 1 , model I 



\z\=-Y^h{J\z\z i+il \) . (8) 



/ . f! 



Since 



1 + Tli ■ Tin 



(9) 



there is thus a ferromagnetic interaction between neigh- 
boring spins that drives an 0(3) ordering transition upon 
increasing J along the JC — axis. The transition re- 
mains in the same universality also for small JC. Since 
the condensate is neutral, the photon remains free and 
is a harmless spectator across this Photon - Molecular 
phase transition. 

Consider now the transition from the Molecular to the 
Higgs phase at large J . This is similar to the one- 
component transition described earlier. Indeed, if we 
take an ordered state of the n's, the remaining U(l) phase 
degrees of freedom of z's are coupled to the gauge field in 
the same manner as the e l6h s in Eq. and the Molec- 
ular to Higgs phase transition is the condensation of this 
charged field. In particular, the transition at J = oo oc- 
curs at the same JC C and is in the inverted 3D XY class; 
the transition is expected to remain in this universality 
also away from the J = oo limit. The Goldstone modes 
of the broken global 0(3) are present both in the Molec- 
ular and Higgs phases and are harmless spectators across 
this transition. 

Having discussed the three phases, consider now the re- 
gion where they meet - the 'fork' in the phase diagram. A 
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b) NCCP1, model I 
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FIG. 1: Phase diagrams of the three models. In the NCCP 
model, * = e ie , while in the NCCP 1 case, * = z = (z T , z{) is 
a two-component field, a) The transition between the Photon 
and Higgs phases is inverted XY and is denoted with open 
symbols, except at JC — oo where it is direct XY (denoted 
with star), b) NCCP 1 model I has a rich phase diagram 
which contains also the Molecular phase where is gapped 
but n orders. Our main interest is on the Photon to Higgs 
phase transition; the points where it is continuous are denoted 
with open circles, while in the vicinity where the three phases 
meet the transition becomes first order and is denoted with 
filled circles, c) The Molecular phase is eliminated in the 
NCCP 1 model II, allowing to focus solely on the Photon - 
Higgs transition. We find that everywhere along the phase 
boundary the transition is continuous (indicated with open 
symbols). 



mean field argument 9 - summarized in Appendix [X] shows 
that in the immediate vicinity of the fork the Photon 
to Higgs and Molecular to Higgs transitions become first 
order, see Fig. [TU The Photon to Molecular transition, 
on the other hand, remains continuous. Focusing on the 
Photon to Higgs transition, the first order line terminates 
at a so-called tricritical point, beyond which the transi- 
tion is continuous. We indeed find similar picture in the 
Monte Carlo simulations of the model, where the transi- 
tion is second-order for the JC > 0.2 points in Fig. [TJ and 
this is the region that we are most interested in. How- 
ever, as we will describe in Sec. IIII El a detailed study of 
the universality is hampered by the fact that this region 
is bordered by two different critical points that are not 
related to the fixed point of our transition: On one hand, 
the system needs to be sufficiently away from the JC = oo 
limit to depart from the 0(4) universality class. On the 
other hand, it also needs to be far from the first-order 
end-point that imposes its own crossovers before the sys- 
tem flows to the fixed point that controls the transition 
we are interested in. Note that the renormalization group 
flows that we have in mind are in a more general space 
than the JC — J parameter space of the model, so the 
above scenario is consistent. 

Because of these complications, the NCCP 1 model I, 
while simple to formulate, requires good understanding 
of all crossovers before interpreting the results. We do 
not have such a complete control to be able to reliably 
extract the critical indices, but we will show represen- 
tative points on the phase boundary where numerics is 
strongly in favor of the second-order character which is 
"more continuous" than in the NCCP case. We will 
further argue for the hypothesis of the tricritical point in 
this model in Sec. Ifff El 

The universality of the Higgs transition that we want 
to study is independent of additional short-range inter- 
actions as long as the transition remains continuous. In 
this respect, to better focus on the scaling properties, we 
also consider the NCCP 1 model II, Eq. ([5]), where we 
have added some antiferromagnetic interaction between 
neighboring spins that hinders the ordering of the n's. 
In fact, with the hindsight of Eq. ([5]), the chosen interac- 
tion precisely compensates the ferromagnetic interaction 
in the original model I along the JC = axis. In par- 
ticular, the model II has no phase transition along this 
axis. We actually find that the Molecular phase is elim- 
inated in this model, and the resulting phase diagram is 
displayed in Fig. [Ib. The transition remains continuous 
all the way to the J = oo limit, and the phase diagram 
in the JC — J plane looks similar to the NCCP case. 

Thus, in the NCCP 1 model II, we avoid the compli- 
cating crossovers caused by the possible tricritical point 
in the vicinity of the Molecular phase. We only need 
to worry about being sufficiently away from the JC = oo 
limit, and this is similar to the NCCP case that we un- 
derstand analytically and can compare numerically. We 
expect that adding generic antiferromagnetic interactions 
between the n's would give a similar improvement in the 
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NCCP 1 model, suppressing the Molecular phase, but we 
have chosen the specific form to avoid possibly introduc- 
ing new phases in the phase diagram. From a Monte 
Carlo study of the model II combined with a separate 
study of the J = oo limit and an analysis of the stability 
away from this limit, we conclude that no other phase 
appears in this model. 

Let us describe in more detail the J ' — oo limit, which 
we formulate as 



yJ=oo 

^NCCP 1 , model 



n [Dz5(\z\ 2 = l) 



exp 



K 



(10) 



Here the primed sum is over 27r-integer-valued conserved 
B-fluxes as in Eq. ([7]). The ensemble also contains CP 1 
fields z residing on the sites and specifying uniquely link 
variables A as follows: 



Ma 



Z - Zn 



(11) 



Note that any 2ir ambiguity in the ^4's can be absorbed 
by the B's. Crudely, we obtain Eq. (fTU|) by focusing on 
the phases of the link variables zjzj while completely 
suppressing the role of their amplitudes, and this is pre- 
cisely achieved when we take the J —> oo limit in the 
model II. 

Comparing with the J = oo limit in the NCCP case, 
the NCCP 1 variant can be also viewed as a loop model, 
but the loops have non-trivial internal structure derived 
from the coupled degrees of freedom Zj. Numerically 
we find that this system has a continuous transition at 
JCc ~ 0.1245, see Sec. IIII C| and the universality of the 
transition is different from the loops with no internal 
structure. We can also carry through an analysis in the 
model II when J is large but finite and find that this 
leads only to short-range modifications that vanish in 
the J — > oo limit (this is similar to how one analyzes the 
NCCP case at large J). We then propose that the whole 
phase boundary in Fig. [TJj is in the same universality as 
the transition found in Monte Carlo at J = oo. 



III. COMPARATIVE STUDY OF THE 
TRANSITIONS 

A. Monte Carlo measurements. Scaling analysis 

We first summarize our numerical measurements. We 
perform classical Monte Carlo simulations using local up- 
dates of the matter and gauge fields, where the latter are 
treated as unconstrained real variables. (We can write 
an explicit finitely defined statistical mechanics model 
that is properly sampled by such simulations, but do 
not belabor this because such details are not used here.) 



In the J — oo models, we also use geometrical worm 
update a 12 ' 13 of the discrete- valued conserved fluxes B. 
Where possible, we use multiple histogram method to 
interpolate between the data points^ the error bars are 
estimated by using blocking method and by running sev- 
eral independent samples. 

As we change the parameters in the system, we mon- 
itor, first of all, thermal properties such as the specific 
heat per site C . This is defined through the variance of 
the action, 



C = 



((S-(S)f 
L d 



(12) 



and can alert about phase transformations in the system 
even when we do not know the underlying physics. At 
a second-order phase transition, C behaves as C ~ \t — 
t c \~ a = \t— t c \~( 2 ~ ud \ where in the last equation we used 
the relation between the specific heat exponent a and the 
correlation length exponent v. The singularity is rounded 
off in a finite system, but manifests itself as a peak in 
the specific heat that evolves characteristically with the 
system size. In a finite volume L d , the peak location 
approaches the bulk critical point as t pea k — t c ~ L^ x l v , 
while the peak height behaves as 



Cm 



C + AL 2/v - d , (second order) 



(13) 



For comparison, at a first order transition, the peak 
height grows with the system size as 

■d 



Cn 



L , (first order), 



(14) 



which would be measured as v r = l/d=l/3ina naive 
analysis in 3D; the peak location converges to the bulk 
critical point as i pca k — t c ~ L~ d . 

Thus, the behavior of the specific heat gives first indi- 
cation about the nature of the transition. In the NCCP 
model, we anticipate a continuous transition with v ~ 
0.67 - the exponent a is then small and negative, and 
the singularity in C is only a cusp. In this case, C is 
less useful for scaling, but a good estimate of v can be 
obtained by studying the third cumulant of the action^ 



C, 



((s-(s)) 3 ) 

L d 



(15) 



which can be viewed as a derivative of C with respect 
to a temperature-like parameter. Near the critical point, 
the characteristic C3 (such as the maximum value or the 
maximum - minimum difference) scales as L 3 /"~ d , which 
is more readily detected. This is useful in the NCCP 1 
model as well, where, previewing the Monte Carlo results, 
the thermal signature of the transition is less singular and 
implies a larger exponent v. 

Let us now consider direct detection of the Higgs tran- 
sition. In the Higgs phase, the gauge field is massive and 
the fluxes do not fluctuate on large scales. On the other 
hand, in the Photon phase, the fluxes proliferate. This 
can be measured by what we call the dual stiffness 



EM*) 



iq-R\2 



<2tt) 2 L 3 



(16) 
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Here the summation is over all plackets oriented perpen- 
dicular to /t, — (V x a) M is the corresponding flux, 
and R is the coordinate of the placket. We consider the 
stiffness at a wavevector q and are interested in the small 
q limit. In an infinite system, pduai is zero in the Higgs 
phase and is finite in the Photon phase. If the fluxes 
are viewed as conserved currents, pduai is essentially the 
superfluid stiffness of this current loop system. 

Recall also that in the NCCP case at J = oo, the 
fluxes b/(2ir) become integer- valued (the corresponding- 
loops are the ANO vortex lines). Pduai is then precisely 
the conventional superfluid stiffness for this loop system 
(which also explains our normalization in Eq. I16[) . The 
loop model is usually studied with no restriction on the 
total flux, and one usually measures the stiffness at q = 0. 
This is not available in our setup: since we are working 
with the gauge fields and periodic boundary conditions, 
the total flux is constrained to be zero. However, we 
can measure the stiffness at a nonzero q ~ 1/L, and in 
the present work we use Pduai and Pduai at the smallest 
9min = (0, 0, 2n/L). In a finite system, this observable is 
less sharp than p q =o- Thus, in the phase with only small 
loops, there is a contribution to p q of order q 2 ~ 1/L 2 , 
while the q = stiffness would vanish exponentially with 
L. Still, we find that this observable works well. The 
NCCP model allows us to develop such experience and 
test our protocols before studying the NCCP 1 case. 

At a second-order transition, we expect the product 
Pduai(Q'min) ■ L to be universal, while near criticality 



Aiuai(<7min; t, L) ■ L = r(f>L 1/l/ ) 



(17) 



where S = t — t c measures deviation from the critical 
point and r(x) is a scaling function. We plot the pduai ■ L 
curves and look at their crossings and thus estimate the 
location of the critical point. We also extract the corre- 
lation length exponent v by applying the above scaling 
form. We try different approaches. One is to consider a 
derivative with respect to a temperature-like parameter 
in the system (such derivatives can be calculated accu- 
rately in the Monte Carlo process). At the critical point 
we expect 



^(Pdual • L) 

dt 1 



L l/u 



(18) 



which can be used to estimate v. This approach is sensi- 
tive to our uncertainty in the location of the critical point 
and is not very robust: e.g., the critical such curve does 
not separate the ones on the two sides of the transition 
but eventually should go above all of them. A different 
procedure is to consider the maximal value of the deriva- 
tive for each L. Yet another approach is to scale the full 
Pduai ■ L curves near criticality. We are using all such 
approaches, but a crude look at the derivatives is simple 
and will already be sufficient for our comparison of the 
NCCP and NCCP 1 models - our main message will be 
that v is larger in the latter case. 

The discussed thermal measures and the dual stiffness 
work in any matter-gauge system. In the NCCP 1 case, 



we also have the field n, Eq. ([6]), which is a local physical 
observable in the matter sector. We measure the magne- 
tization associated with the 0(3) ordering, 



(19) 



and estimate the critical point from the crossings of the 
corresponding Binder cumulant ratio 



I M I 



(\M\ 



2\2 



(20) 



Similarly to Eq. (fl7|) , we expect the following scaling 
form near the critical point: 



g{t,L)=g{8L 1 '») 



(21) 



which we analyze in the same spirit as the pduai ■ L. 

We characterize the 0(3) ordering by the magnetiza- 
tion exponent (3: m ~ 5@. This gives the following finite- 
size scaling ansatz 



m(t,L) = L-PlvfiSL 1 '") 



(22) 



In particular, right at the critical point we expect 

mU-I-"/^^)/ 8 , (23) 

where the last equation defines exponent r\ in 3D. By 
plotting to vs L on a log- log plot, the critical curve is 
a straight line that separates the curves on the ordered 
and disordered sides. This analysis independently checks 
the location of the critical point and gives an estimate of 

Finally, we also observe the transition by monitoring 
the hclicity modulus T for twisting the 0(3) vector field 
ft. This can be calculated in the standard way during 
the same Monte Carlo process, and we give detailed ex- 
pressions for the two NCCP 1 models in Appendix [B] 
The helicity modulus is also a measure of the ordering 
in the 0(3) matter sector, but has a somewhat different 
character from the magnetization and is more like a su- 
perfluid stiffness in a boson system. The scaling of the 
hclicity modulus at a second-order transition is similar 
to the dual stiffness described earlier. In particular, we 
expect universal crossings of the T • L curves, which is 
yet another way of locating the transition in our finite- 
size study. These measurements can also be compared 
with other studies of such matter-gauge systems^ and 
also with studies of quantum Hamiltonians whose criti- 
cal behavior may be described by such field thcories4££ 



B. Representative study of the transitions in the 
three models at K=0.4 

Here we present Monte Carlo results at JC = 0.4 in the 
three models. This is sufficiently away from the K, = oo 
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limit - for example, the critical coupling J C (JC = 0.4) 
is at least 50% larger than J C {JC = oo) in each case: 
Jc(0.4) = 0.682 vs J c (oo) = 0.4541 in the NCCP model; 
J c (0A) = 1.412 vs Jc(oo) = 0.9358 in the NCCP 1 model 
I; and J~ c (0.4) = 1.68 vs J c (oo) = 0.98 in the NCCP 1 
model II. Furthermore, the transition is clearly observed 
in the behavior of the gauge field as detected by the dual 
stiffness - we see a convergence of the crossings of pduai ■ L 
for our system sizes. In the NCCP 1 models, we also see 
convergence of the Binder ratio crossings and the T • 
L crossings, which detect the transition in the matter 
sector. Recall that the gauge sector is completely absent 
in the JC = oo limit, and indeed at large JC the "critical 
points" defined by the "gauge" and "matter" measures 
start far apart for our system sizes, but are already close 
at JC = 0.4. We expect comparable finite size effects in 
the three systems as far as the development of the true 
Higgs criticality is concerned. 

We will discuss the observed trends along the Photon - 
Higgs phase boundary in Secs. lIIID|IIII El The represen- 
tative point JC = 0.4 was chosen to be sufficiently away 
from the first-order region in the NCCP 1 model I, with 
the hope that it is not affected strongly by crossovers 
near the tricritical point. We argue here that the transi- 
tion is continuous in all three models at JC = 0.4. In the 
NCCP case, this is expected theoretically and is sup- 
ported numerically. In all measured respects, the tran- 
sition in either of the NCCP 1 models at JC — 0.4 looks 
more continuous than NCCP . 

Figure [2] shows the specific heat per site in the three 
models. In each case, we see a peak whose position moves 
towards the estimated critical point and whose height in- 
creases gradually. The slow evolution of the position and 
the weak (at most sublinear) growth of C max with the 
system size are indicative of the continuous nature of the 
transition and rather large correlation length exponent 
v in each case. Fitting C max to Eq. (flU)) . we estimate 
v w 0.65 — 0.7 in the NCCP case, which is close to the 
accepted value in the (inverted) 3D XY. We find a some- 
what larger v w 0.70 in the NCCP 1 model I case. In 
the NCCP 1 model II, the growth of the peak height with 
the system size is very weak and we do not attempt to 
analyze C max . We also measure the third cumulant C3, 
Eq. (|15p , and analyzing this data gives estimates of v con- 
sistent with the above. In the three systems at JC = 0.4, 
we do not see any tendency towards first order behav- 
ior such as Eq. (|14[> . Furthermore, energy histograms for 
system sizes up to L = 36 show one mode and no sign 
of any structure developing. From all such measures, the 
transition in the presented NCCP 1 systems has weaker 
thermal singularities than in the NCCP case. 

Let us describe how we locate the transitions more ac- 
curately. Figure [3] shows the crossings of pduai(<7min) ■ L, 
which approach the critical point from the opposite di- 
rection compared with the C max positions. The crossings 
converge quickly for our system sizes and thus allow an 
accurate determination of J c . 

In the NCCP 1 systems, we also observe the transition 
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FIG. 2: (color online). Specific heat C in the three mod- 
els at JC — 0.4 varying J near the transition. System sizes 
are L — 12 — 32. In the top and middle panels, the lines are 
obtained by multiple histogram method, which provides a sys- 
tematic interpolation between data points; the corresponding 
peak values C max are plotted in the insets and show sublin- 
ear growth with L. In the bottom panel, the lines are only 
guides to the eye connecting nearby data points; in this case, 
the peak value changes too little to be extracted reliably. In 
each panel, the vertical arrow shows the bulk critical point J c 
estimated from Figs. [3] and [4] The specific heat peak moves 
towards this point gradually and the developing singularity is 
weak, which is a direct evidence of the second-order nature of 
the transition. 
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directly in the matter sector by measuring the magne- 
tization to, Eq. (|19[) . The corresponding Binder ratios 
are shown in Fig. [4] for the models I and II. Comparing 
with the appropriate panels in Fig. [3J we see that the 
Binder crossings move in the opposite direction to the 
/9duai(<?min) ' L crossings, and the two measures effectively 
bound the location of the transition to a narrow region. 

Finally, Fig. [5] shows the crossings of T • L in the two 
NCCP 1 models, providing an alternative detection of the 
transition in the matter sector. The crossings converge 
towards the bulk critical point similarly to the Binder 
crossings and also help to bound the transition region. 

Having estimated the critical points, let us now discuss 
scaling in the three models. In the panels in Figs. [3J |U 
and O we also show the corresponding derivatives eval- 
uated at points bounding the transition^ vs the system 
size L. In each case, we perform simple linear fits cor- 
responding to Eq. I|18p on the log-log plot and list the 
extracted values of v. This procedure is not very reli- 
able and does not necessarily bound the exponent, but it 
gives a crude picture of the results we obtain doing more 
elaborate scalings. 

In the NCCP case, such scalings of Pduai(<Zmin) ■ L 
suggest v between 0.5 and 0.55, which we interpret as 
some effective value on our length scales. This is signif- 
icantly different from the expected v ~ 0.67, but when 
we look closer, we do see the corresponding plots curving 
in the direction of such larger v. Thus, naive scaling of 
Aiuai(<7min) -L in the one-component matter-gauge system 
at this JC suffers from strong finite size effects. We could 
in principle improve our estimates of v by performing 
more complex fits including corrections to scaling, but in 
the absence of good control over these, we do not venture 
in this direction. Instead, in the spirit of the compara- 
tive study, we should be prepared for similar effects also 
in the two-component system. 

In the NCCP 1 model I, such scalings of /9duai(9min) 1 L 
show effective v between 0.55 and 0.65. This is larger 
than the effective v in the above NCCP case, suggest- 
ing that the true v in the two-component system is larger 
than in the one-component system. From the scaling 
of the Binder ratio, of the T • L, and of other observ- 
ables in the matter sector, we estimate v between 0.65 
and 0.7 (somewhat smaller than the crude values in the 
top panel in Fig. 2]), which further supports our proposal 
^nccp 1 > ^nccp ■ We warn, however, that our scalings 
in the model I at fC = 0.4 may already be affected by the 
proximity of the Molecular phase - see the discussion of 
the trends along the Photon - Higgs phase boundary in 

sec. nnii 

Turning to the NCCP 1 model II, the scalings of 
Pduai(<7min) • L consistently suggest v around 0.7 (some- 
what larger than the crude estimates in the bottom panel 
in Fig. [3J), while the scalings of the Binder ratio suggest 
v between 0.7 and 0.75. This system does not have any 
Molecular phase nearby and, being sufficiently away from 
the K. = oo limit, is closer to the true two-component 
Higgs criticality than the model I. The T ■ L data in the 
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FIG. 3: (color online). pduai(<?min) ■ L in the three models 
at K, = 0.4. The crossings converge quickly towards J c = 
0.682 in the NCCP model (top panel), J c = 1.412 in the 
NCCP 1 model I (middle panel), and J c = 1.68 in the NCCP 1 
model II (bottom panel). In the last two cases, we can check 
the convergence against the Binder ratio crossings, Fig. [4] 
which move in the opposite direction. In each panel, the inset 
summarizes naive estimates of the correlation length exponent 
v applying Eq. (|18|l to the derivative data at points bordering 
the transition - see text for details. 
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FIG. 4: (color online). Binder ratios in the NCCP 1 models 
I and II at K, = 0.4. The crossings bound the transition 
point from the opposite side compared with the Pduai(Qmin) ■ L 
crossings, Fig. [3] In the insets, we show naive estimates of v 
fitting the derivative of the Binder ratio to the form ~ L 1 /" 
at points bordering the transition. 
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FIG. 5: (color online). T • L in the NCCP 1 models I and 
II at K, = 0.4. The crossings move in the same direction as 
the Binder ratio crossings, Fig. [4] and bound the transition 
point from the opposite side compared with the Pduai(gmm) ■ L 
crossings, Fig. [3] The inset in the model I case shows scaling 
analysis of the derivative similar to the analysis in Fig. [3] (no 
such derivative data was measured in the model II). 



model II is rather noisy, as one see from the correspond- 
ing panel in Fig. This model is harder to simulate 
than model I: e.g., we cannot use local heat bath updates 
of the matter fields and we cannot use multi-histogram 
reweighting; there is also significant complexity in the 
expression for T, Eq. (|B1[) . Still, we see that T • L shows 
reasonable crossings and is bounded by value of order 1 
at criticality with no significant drifts. Here, it is helpful 
to have the bounds on the critical J c from the more accu- 
rate Pduai • L and Binder crossings, Figs. [3] and 01 so there 
are no uncertainities where the T -L crossings might flow. 
We will describe overall trends along the phase boundary 
in the model II in Sec. ITTTDl 

Finally, from the magnetization scaling, we also es- 
timate the exponent rj in the two NCCP 1 models at 
JC = 0.4, finding it to be around 0.4 — 0.5 in both cases 
on our system sizes. 



C. NCCP model and NCCP 1 model II at J -> oo 

Before describing overall trends along the Photon - 
Higgs phase boundary in the three models, we present 
Monte Carlo results in the two J — oo models, Eqs. ([7]) 
and (jlOp . In our simulations, we also use geometri- 
cal worm update s 12 i 13 of the discrete-valued conserved 
fluxes B. This eliminates the critical slowing down in 
the NCCP case and significantly improves statistics in 
the NCCP 1 case. In the latter, we still use local updates 
of the matter fields, which is probably why the critical 
slowing down persists. 

As we have already described, the NCCP model at 
J = oo is equivalent to the well-studied loop model with 
short-range interactions. We use it to develop experience 
with the Pduai(<7min) • L to characterize the transition, 
and we find that this measure works well. The study 
also allows us to connect with the points measured along 
the phase boundary in the NCCP model at finite J , 
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Fig. [T^,. gaining experience with finite size effects in such 
matter-gauge simulations - we will use this when dis- 
cussing overall trends in Sec. IIIIDI Furthermore, we can 
now compare the two J = oo models Eqs. ([TJ and (fTU|) 
directly. 

In Fig. [SI we show the specific heat C in the two mod- 
els. In the NCCP case, there is a clear peak which moves 
towards the bulk critical point and sharpens slowly with 
increasing system size. Analyzing the thermal singularity 
as in Sec. IIII Bl and as described in Eqs. fT2"|) - (fl"5")) . we ob- 
tain an estimate of the correlation length v = 0.66 — 0.67, 
which is in agreement with the accepted value v = 0.67. 
Turning to the NCCP 1 case, the specific heat signature 
of the transition is weaker and appears as a feature that 
develops on top of a sloping down background. The evo- 
lution of the feature with L is more slow than in the 
NCCP case, implying a larger exponent v. From the 
study of the third cumulant of the action, we estimate 
v~ 0.75 -0.8. 

Let us now characterize the transitions using the dual 
stiffness shown in Fig. [71 In either case, we see a rapid 
convergence of the Pduai(<Zmin) • L crossings towards the 
putative bulk critical point. In the NCCP model, the 
crossing for the largest two systems L — 28 and 32 
occurs around JC = 0.07609, very close to the true 
JC C = 0.07606, and the crossing location is still moving 
slightly in the expected direction. In the inset, we show 
the derivative rf(paual • L)/dJC as a function of L at two 
points JC that bound the critical point; direct fits give 
v&t (K = 0.0760) = 0.62 and u^QC = 0.0761) = 0.73, 
which in turn bound the expected v. As discussed in 
Sees. IIII Al and [Til Bl the inset should be viewed only as 
a crude illustration of the scaling analysis, but it is al- 
ready sufficient for our comparative study. 

Proceeding with similar analysis in the NCCP 1 case 
shown in the bottom panel in Fig. [71 the critical point 
is between JC = 0.1240 and 0.1245. From the analysis in 
the inset, we estimate v between ^fit(^C = 0.1240) = 0.68 
and ffit(/C = 0.1245) = 0.80. More elaborate scaling 
procedures suggest that v is between 0.7 and 0.75. The 
correlation length exponent is thus detectably larger than 
in the NCCP case. 

In the NCCP 1 case, we also consider the Binder ratio, 
Fig. [SI The crossings of g bound the critical point from 
the opposite side compared with the pduai ■ L crossings. 
Performing finite size scaling analysis of the Binder ratio 
curves gives estimates of v between 0.7 and 0.8, which is 
consistent with the presented other estimates. 

Finally, in Fig. [51 we show finite size scaling analy- 
sis of the magnetization. We plot m(JC, L) as a func- 
tion of L at fixed JC, expecting Eq. ([23)) at the criti- 
cal point. In the figure, there is a clear separation be- 
tween the magnetically ordered and disordered sides, and 
from the fits at the test points we estimate r\ between 
Vm (JC = 0.1240) = 0.29 and T] m (JC = 0.1245) = 0.18. 
Since the Binder ratios, which detect the criticality in the 
matter sector, still cross near JC = 0.1240 for our largest 
systems, the correspond value of rj is likely closer to the 
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FIG. 6: (color online). Top panel: Specific heat in the NCCP 
model at J — oo, Eq. (fTJ. This is the same as a loop model 
with steric interactions in an ensemble with no windings. Bot- 
tom panel: Specific heat in the NCCP 1 model II at J — oo, 
Eq. (|10p . In each panel, the connecting lines are obtained by 
multiple histogram method; the vertical arrow shows the bulk 
critical point JC C estimated from Figs. [7J and [8] The NCCP 1 
transition sharpens more slowly with L than NCCP and is 
characterized by a larger exponent v. 



true exponent. When we try to scale the full magneti- 
zation curves near criticality, we get a somewhat wider 
range of r\ = 0.2 - 0.4. 

To summarize, the transition in the NCCP 1 model II 
at J = oo has a weaker thermal singularity than in 
the NCCP model; it is a second-order transition with 
^nccp 1 > ^nccp - This conclusion is also supported by 
the scaling analysis of the pduai ■ L and of the Binder ra- 
tios, which give 

v = 0.7 -0.75, (NCCP 1 model II, J = oo). (24) 

The magnetization order is characterized by the exponent 

7? = 0.2 -0.4, (NCCP 1 model II, J = oo). (25) 
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FIG. 7: (color online). Top panel: Pduai(9min) ■ L in the 
NCCP model at J — oo. The critical point is known ac- 
curately to be IC C = 0.07606 [y — 0.67), and the crossings 
converge rapidly to this location. Inset shows the derivative 
of pduai • L with respect to K, as a function of L at two test 
points that bound the critical region; the lines are direct fits 
to the form ~ L 1 ^ , and the corresponding estimates of v are 
also listed. Bottom panel: Similar analysis in the NCCP 1 
model II at J = oo; our best estimate for JC C is K, c ~ 0.1245. 



1. T ■ L in the NCCP 1 model II at large J = 16 

We have not been able to measure the helicity modulus 
T directly in the NCCP 1 model II at J = oo: One ex- 
pression for the T that we derived in this limit contained 
singular integrals over the z variables and was ill suited 
for the Monte Carlo sampling integration. Instead, wc 
performed a study at large but finite J = 16 and mea- 
sured the expression Eq. (|B1|) . This is our closest con- 
nection with the J = oo limit tracking the behavior of 
T-L along the phase boundary in Fig.[T}:. The horizontal 
J = 16 cut is very far from the decoupled K, = oo limit 
and also far from the vertical K. = 0.4 cut presented in 
Sec. IIIIBI Nevertheless, the measured T • L crossings on 
the two cuts are rather similar, illustrating that there is 
no significant flow of the critical T • L values along the 
phase boundary. 
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FIG. 8: (color online). Binder ratio g in the NCCP 1 model II 
at J — oo. Comparing with the corresponding pduai ■ L plot 
in Fig. the g crossings converge to the same critical point 
bounding it from the opposite side. Inset shows estimates of 
v from the analysis of the derivative of g, in the same spirit 
as in Fig. [7J 
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FIG. 9: (color online). Finite size scaling analysis of the 
magnetization. We expect m ~ j / ~( 1 +'')/ 2 a t criticality and 
perform fits for the exponent r\ at two test point. On the log- 
log plot, the ordered and disordered sides curve up and down 
respectively, confirming the estimates of the critical region 
from Figs. [7] and HI 



The J = 16 study is summarized in Fig. |TD] The 
/?duai(<7min) • L crossings and the Binder ratio crossings 
behave similarly to the J = oo model, Figs. [7] and [S] 
These measures, which have smooth limits for large J ^ 
have clean statistics and allow us to accurately constrain 
the critical point to lie between JC = 0.132 and 0.1325. 
On the other hand, the T • L data is much more noisy: 
Individual terms present in Eq. (|B 1 1) can have large typ- 
ical values and large cancellations are to occur upon the 
full integration. Still, we see that the T • L values near 
the criticality have not changed significantly from the 
K, = 0.4 cut in Fig. despite the much larger value of J 
here. 
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FIG. 10: (color online) Universal crossings in the NCCP 1 
model II at large J = 16. This is our closest approach towards 
Cf = oo, and the focus here is on the boundedness of the 
T • L at criticality. In the top and middle panels, pduai(gmm) ■ 
L crossings move towards smaller JC values and the Binder 
ratio crossings move towards larger JC values, localizing the 
transition between JC = 0.132 and 0.1325. The T • L data 
in the bottom panel is noisy, but it is bounded by a value 
close to 1 and is similar to the T • L crossings at other cuts 
across the phase boundary in Fig.[TJj (compare, e.g., with the 
vertical cut at JC = 0.4 shown in Fig.[SJ). In the top and middle 
panels, where the data quality is good, the interpolating lines 
are obtained by the multi-histogram reweighting technique, 
while in the bottom panel the lines connect the individual 
data points; the error bars are estimated by averaging over 
different Monte Carlo runs. 



D. Qualitative summary of the full Photon - Higgs 
phase boundary in the NCCP model and NCCP 1 
model II 



We now summarize Monte Carlo results along the Pho- 
ton - Higgs phase boundary in each of the three models. 
The phase diagrams were mapped out using the same 
criticality locators as in the preceding sections. Small 
systems up to L = 16 — 20 were sufficient to get good 
estimates of the transition points marked in Fig. [TJ while 
several points were studied up to L = 24 — 32 and with 
detail comparable to Sees. IIII BIHI Cl 

Specifically, in the NCCP model and in the NCCP 1 
model II, we studied a vertical cut at JC = 0.25 and a hor- 
izontal cut at J = 4, cf. Figs. [T^,,c. Combined with the 
study of the J = oo limit, we thus observed the full evo- 
lution of our measurements along the phase boundary in 
the two models. Some details will be given shortly, but 
we can already say that in either model the criticality 
at different places along the phase boundary is qualita- 
tively similar to the presented JC = 0.4 case and connects 
smoothly to the J = oo limit. From all studies, the 
J = oo models are probably closest to the desired Higgs 
universalities, being least affected by other fixed points. 
In particular, these models provide our best estimates of 
the critical indices. We also studied in detail the phase 
boundary in the NCCP 1 model I, where the situation is 
more complicated - we will return to this model in the 
next section. 

First, we give some details in the NCCP case. The 
transitions are detected using the pduai ■ L crossings and 
look qualitatively similar to the JC = 0.4 case presented 
in Sec. iHlBl and the J = oo case in Sec. InTCl The 
locations of the crossings converge quickly to the critical 
point. On the other hand, the specific heat peak evolves 
slowly with the system size and moves towards the bulk 
critical point as in the top panels of Figs. [2] and O 

One observation that we make is that the values of 
Pduai • L at the crossings, while apparently converged, are 
different at different locations along the phase boundary 
- compare, for example, the top panels of Figs. [3] and 0- 
which is at first disconcerting. However, the observed 
values at the crossings approach those at J = oo as 
we move towards this point along the phase boundary. 
We believe that in the thermodynamic limit, /?duai ■ L is 
universal at criticality. Indeed, we can show that our 
dual stiffness in the one-component matter-gauge sys- 
tem is equal to a sum of a stiffness p of a mathemati- 
cally equivalent discrete-valued current loop model with 
short-range interactions and of a positive contribution of 
order \j(JJL ). Then, assuming universal p-L at critical- 
ity in discrete loop models, the pduai • L crossings in the 
NCCP model must converge to the same value but with 
a positive finite-size correction of order 1/(J7X), which is 
indeed comparable to the magnitude of the difference be- 
tween the pduai ■ L values in Figs. [3J and [7] We note also 
that in the NCCP model at K = 0.4, the initial drift of 
the pduai • L value at the crossings upon increasing our 
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system sizes is actually away from the J = oo value, but 
we expect that this drift will eventually turn back. While 
we do not have full understanding of all such finite-size 
effects in our systems, the smooth evolution along the 
phase boundary and our analytical understanding make 
us confident in the proposed inverted 3D XY universality 
in the NCCP model. It is also useful to note that such 
corrections will have much smaller effect on the rapid 
convergence of the locations of the crossings to the bulk 
critical point, which is indeed what we observe. 

Concluding with the NCCP case, the scaling analysis 
of the dual stiffness as in Sec. IIII Bl yields effective v that 
starts around 0.5-0.55 for JC = 0.4 and moves towards the 
expected value v = 0.67 as we approach J = oo, adding 
to our confidence that we understand this system well. 

Let us now provide similar details for the NCCP 1 
model II, whose phase diagram has the same topology 
as in the NCCP case. Here again we see no qualitative 
difference in the observed criticality as we move along the 
phase boundary from the JC = 0.4 point to the J = oo. 
For example, the specific heat peaks evolve slowly to the 
bulk critical points as in Figs. [2] and [6l The peak heights 
C m ax in the vertical cut at JC — 0.25 and in the horizontal 
cut at J = 4.0 are essentially size- independent, i.e., they 
look even less sharp compared to the JC = 0.4 case. Us- 
ing pduai • L to detect the criticality, the picture that we 
observe is also very similar to the NCCP model. In par- 
ticular, since the two systems are simulated under similar 
finite size conditions, it is likely that similar effects are 
responsible for the apparent difference in the pduai ■ L val- 
ues at the crossings as we move along the phase boundary. 
In this respect, the NCCP model, which we otherwise 
understand well, provides a useful reference on such ef- 
fects that we still don't have much experience with. This 
gives us confidence when interpreting the NCCP 1 case, 
which is our main focus. 

In the NCCP 1 model, we also observe the transition 
using the ordering of ft, which monitors the matter sector. 
In a broad view of the phase diagram, Fig. [TJ;, for large 
JC the Binder ratio crossings and the pduai ■ L crossings 
start far apart for our system sizes, but as we move along 
the phase boundary towards smaller JC, the separation 
in the apparent criticality in the two sectors decreases 
significantly. Recall that the gauge sector is absent in 
the JC = oo limit. The convergence of such locators of 
the criticality in both sectors shows that the system is 
sufficiently away from the JC = oo limit and we are indeed 
probing the true Higgs universality. 

At points between JC = 0.4 and J — oo, scaling anal- 
yses of both the pduai ■ L and the Binder ratio give us 
consistently v in the range 0.7 to 0.75, supporting the 
conclusion that we arc dealing with the same continuous 
transition for all points on the phase boundary in this 
model. Finally, our study of the helicity modulus also 
agrees well with the pduai ■ L and Binder ratio measure- 
ments. Thus, the T • L crossings agree with the other lo- 
cators of the criticality. Furthermore, everywhere on the 
phase boundary, wc find roughly similar values T-Lsl 
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FIG. 11: (color online). Peak height of the specific heat in 
the NCCP 1 model I measured at different locations along the 
Photon - Higgs phase boundary, cf. Fig. [lb. The data at 
JC = 0.4 is the same as in the corresponding inset in Fig. [2] We 
remove the overall scale of the action by plotting C mBX _/J^, 
where J c is the critical value for a given JC. Similar plot in 
the NCCP 1 model II (not shown) has C m ax that depends very 
weakly on L for all points along the phase boundary in Fig.QJ. 

at criticality, including the cut at J = 16 presented in 
Sec. IIII C 11 which is the largest J for which the T was 
measured. There is no significant flow of the critical T ■ L 
values along the phase boundary in the NCCP 1 model 
II, consistent with the same continuous transition every- 
where as implied by all other measures. 



E. Qualitative summary of the full Photon - Higgs 
phase boundary in the NCCP 1 model I 

Let us now discuss the NCCP 1 model I. Fig.llllsumma- 
rizes our measurements of the thermal properties along 
the Photon - Higgs boundary, cf. Fig.[lJ>. For JC > 0.25, 
the specific heat peak height C max grows weakly with 
L - at most sublinearly for our system sizes. For JC = 
0.18 — 0.20, the growth is linear, and becomes faster and 
dramatic for smaller JC. Note that in Fig. [IT] we have 
divided C max by the corresponding thus removing 
crudely the overall scale, so it is reasonable to compare 
the magnitudes among the different JC data. The pduai • L 
crossings, Binder cumulant crossings, and T • L cross- 
ings are already well converged for these system sizes 
(cf. the JC = 0.4 data in Sec. IIII Bp . i.e., these systems 
are far from the JC = oo limit and represent fully cou- 
pled matter-gauge criticalitics. Our interpretation of the 
thermal signatures is that the transition is second-order 
for JC > 0.2 and is first order for JC < 0.2, while JC « 0.2 
is a candidate for the tricritical point separating the two 
behaviors. (It is worth repeating here that a similar C max 
plot in the model II, where there is no Molecular phase, 
has essentially size-independent C max for all points along 
the phase boundary.) 
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The transition at JC = 0.3 appears more sharp com- 
pared to the JC = 0.4 point presented in Sec. IIII Bl 
and various broad-brush scaling analyses give a smaller 
v c fi{K, = 0.3, L < 32) as 0.6. The transition appears 
still more sharp at IC = 0.2 with v e s sa 0.5. Also, the 
effective exponent r\ for the magnetization order param- 
eter decreases towards zero as we approach JC = 0.2. For 
K, = 0.25,0.3, we observe that the extracted v and r] 
tend to increase slightly when the scaling range around 
the critical point is made more narrow. We interpret the 
apparent exponent drift along the phase boundary to be 
due to crossovers near the tentative tricritical point near 
IC = 0.2. (To repeat here, we see no such exponent drift 
in the model II). 

Certainly, this drift is an important concern that we 
have been unable to fully address in the model I with our 
system sizes. These difficulties are the main reason why 
we also considered the NCCP 1 model II, where the first 
order region is absent and the systematics of the scalings 
is improved, as we have already described. This experi- 
ence with the model II gives us better grounds for making 
the above interpretation in the model I. It would be use- 
ful to confront this with direct simulations in the model I 
on larger system sizes than in our work. Larger systems 
should allow one to better filter out the crossovers and 
focus on the true critical behavior, e.g., confronting our 
expectation of the universal exponents along the contin- 
uous region of the phase boundary. 

To give more details behind the proposed scenario in 
the NCCP 1 model I, we also present the behavior of the 
Binder ratio crossings and T • L crossings along the Pho- 
ton - Higgs phase boundary. These are vital locators of 
the criticality and are expected to converge to universal 
values in the second-order case. Going back to Figs. 0] 
and [5] at JC = 0.4, we see some weak evolution of the 
values at the crossings. We find that these values further 
drift as we move along the phase boundary. Similarly to 
the drifts in the apparent exponents, this is an important 
concern. In this regard, it is helpful to have a picture of 
all data along the phase boundary. We can then try to 
see possible attracting points of such drifts of the cross- 
ings and also see where the presented detailed study at 
IC = 0.4 of Sec. IIII Bl stands in the overall picture. 

The top panel in Fig. [T2l summarizes all collected data 
for the Binder ratio crossings. For each JC, we consider 
a vertical J scan across the phase boundary in Fig. [TJa 
and analyze the Binder ratios as done in Fig. |4j For 
each pair of system sizes L\ and L2, we find the corre- 
sponding crossing of the curves and organize the results 
by showing {y/L 1 L 2 , g cmss {L\,L2\JC)}. Admittedly, the 
plot Fig. [12] is noisy, since the crossings are fairly sensi- 
tive to the statistical errors. For some systems we show 
results obtained from several independent data samples, 
so the spread of the symbols gives some idea about our 
uncertainties. Generally, we do not trust the values when 
the crossings become non-systematic with L, which hap- 
pens for our largest systems. From the available data, 
we see that for JC > 0.2 the Binder ratio values evolve 
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FIG. 12: (color online). Top panel: Binder ratio values at the 
corresponding crossings for all pairs of system sizes L\ and 
L2. The data is collected by analyzing J scans such as Fig. [4] 
at a number of IC locations along the Photon - Higgs phase 
boundary from Fig. [lb. For an illustration, the pairs of system 
sizes are organized by \JL\Li on the horizontal axis. Bottom 
panel: T • L values at the corresponding crossings obtained by 
analyzing the helicity modulus plots such as Fig. [5] Together 
with Fig. [T3] and all other data, our interpretation is that the 
transition is continuous for IC > 0.2 and becomes first order 
for JC < 0.2. 



strongly initially but then appear to saturate to a value 
g m 1.2 for the largest systems. On the other hand, for 
JC < 0.2 the Binder ratios at the crossings flow to smaller 
values. 

Note that in the scenario with the tricritical point, we 
expect that the Binder ratios go to different nontrivial 
values in each of the three cases - at the second-order line, 
at the tricritical point, and at the first-order transitions. 
In this respect, the fact that the initial evolution in all 
cases is in the same direction makes the interpretation 
of the Binder crossings data more difficult. Certainly, 
improving the statistics of the data in Fig. [T^] and going 
to larger sizes could address such concerns and confront 
our interpretation, which we nevertheless make given all 
available data. 

A more clear separation between the different regimes 
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is provided by similarly collected data for the T • L cross- 
ings shown in the bottom panel in Fig. 1121 In the second- 
order case, we expect the crossings to go to a universal 
value, and we see that for IC > 0.2 these indeed appear to 
accumulate to a value T • L w 1 for the largest systems. 
On the other hand, in the first-order regime, we expect 
the T • L crossing values to increase linearly with the sys- 
tem sizes, and we indeed see such growth for IC < 0.2. 
We are fairly confident that in the proposed second-order 
region the values of T • L do not grow without limit at 
the transition, since we can bound the possible drift of 
the T • L crossing locations by the pduai ■ L crossing lo- 
cations that converge to the transition from the opposite 
side, cf. Figs. [3j [U and O (We repeat here that the 
T • L crossings in the NCCP 1 model II were also found 
to be bounded and of order 1 at several different places 
along the phase boundary; the situation in the proposed 
second-order regime in the NCCP 1 model I is consistent 
with this in the NCCP 1 model II.) 

In Fig. ll3i we present a somewhat different perspective 
on the collected Binder ratio, T-L, and pduai ■ L data. For 
each /C, we define an operational critical point J C {V) for 
a given size L as the point where the Binder ratio takes 
value g = 1.15. We then measure T • L and pdual ■ L 
at this point and plot the results as a function of L. 
The specific value g = 1.15 lies between the apparent 
crossing values in the IC > 0.2 and IC < 0.2 regimes 
and is selected so that J C (L) approaches the true critical 
point systematically from the large J values when IC > 
0.2. Any g between the limits on the two sides of the 
transition can be used for such a definition of the finite- 
size critical point, and we have checked that the results 
do not change qualitatively if we fix the Binder ratio at 
different values such as g — 1.3, 1.2, 1.1 instead. 

From the top panel in Fig. [T31 we see that thus mea- 
sured T • L values hardly change with the system size for 
IC > 0.2 and all are in a narrow range around 1.2 — 1.3. 
There may be slight systematic drift, and the ultimate 
fate of this is hard to judge given the uncertainties of 
our measurements and limited system sizes. However, 
we note that the relative variations of the T • L values 
are small: Had we chosen the vertical axis in the figure 
to start at zero, which is natural given that there is no 
other scale for the expected universal values of T • L, the 
curves for JC > 0.2 would look essentially flat. More- 
over, the values are fairly close in a whole range range 
IC = 0.25 — 0.6, and are also similar to the values we find 
at this Binder ratio in the NCCP 1 model II along the 
whole phase boundary (e.g., near points IC — 0.4,0.25 
and J = 4.0, 16.0 in Fig. [lb). This is to be contrasted 
with the behavior for IC < 0.2, where we see a strong 
growth of the critical T • L with the system size start- 
ing from the small L, which is interpreted as first-order 
behavior. 

Similar considerations apply to the pduai • L measured 
at fixed g = 1.15 and shown in the bottom panel in 
Fig. [T3J For IC < 0.2, the pduai ■ L values grow strongly 
with the system size, while for IC > 0.2 they grow very 
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FIG. 13: (color online). Different view of the Binder ratio, 
T • L, and /5d ua i ■ L data. For each system size, we define the 
finite-size critical point to be where the Binder ratio takes 
value g = 1.15 and then measure the T • L and Pduai ■ L at 
this point and plot the results as a function of L. Our inter- 
pretation of the data is that the transition is continuous for 
IC > 0.2 and becomes first order for IC < 0.2. 



weakly and tend to saturate. We should mention that 
we do not see such clear separation between the two 
regimes by looking at the pairwise crossings of pduai • L 
(not shown): The values at the crossings are essentially 
size-independent for all IC for our sizes, i.e., we do not 
see the first order region in this measure yet. We also 
want to point out that from the study of the model II at 
J = oo, Fig. [3 we expect the critical value of pduai ■ L to 
be eventually smaller than what we see at the crossings 
in the model I. That is, there are strong finite size correc- 
tions to this measure, but such tendencies are similar to 
what we have observed in the NCCP model and NCCP 1 
model II and discussed earlier. 

We have also analyzed the evolution of the extracted 
T] exponent in the two NCCP 1 models along the phase 
boundary. In the model II, we see a convergence with 
the system size towards the value r\ ss 0.3 estimated in 
the J — > oo limit, cf. Eq. (J2H). On the other hand, the 
effective r\ is still larger around 0.4-0.5 in the model I 
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at JC = 0.6, 0.4 for our system sizes but is decreasing. 
The initial large effective 77 is expected, given that this 
'Molecular' order parameter has r\ sa 1.37 at the 0(4) 
pointii where the matter and gauge fields are decoupled. 
This also explains why the original small system study of 
the NCCP 1 model I at JC = 0.6 in Ref. H obtained large 
77, which we now revise. 

To conclude, despite some uncertainties with the weak 
drifts as we move along the Photon - Higgs phase bound- 
ary in the NCCP 1 model I, we believe that there is a 
wide regime where the transition is continuous. While 
the discussion of the full phase boundary is sketchy (in 
particular, we have not tried to accurately characterize 
the tricritical and first-order regimes), it provides bet- 
ter understanding of possible issues and concerns in our 
numerical explorations of the Higgs transition, and, in 
particular, gives a better idea where the detailed study 
at K = 0.4 of Sec. IIII Bl stands in the overall picture. 
We conclude that JC = 0.4 is well inside the continuous 
regime; recall that it is also far from the JC = 00 limit, so 
we have strong reasons to claim that it is representative 
of the true continuous Higgs criticality. 



IV. SUMMARY AND CONCLUSIONS 

To summarize, we perform a comparative Monte Carlo 
study of the Higgs transition in the one-component and 
two-component lattice superconductors and argue for the 
generic second-order nature in both cases. The one- 
component model (called NCCP here), which is well 
understood analytically and numerically, serves as a ref- 
erence system for our matter-gauge simulations. 

We first examine in detail the two-component model 
introduced in Ref. 2 and called NCCP 1 model I here. The 
phase diagram is shown in Fig. [JJd, and we argue that 
there is a range of parameters, JC > 0.2, where the tran- 
sition is second-order. The main evidence is the weakness 
of the thermal singularities, but other measures are also 
in favor of the second-order nature. We present an ex- 
ample of such an analysis at JC = 0.4 and argue that the 
correlation length exponent is larger than in the NCCP 
case. 

The NCCP 1 model I also has the Molecular phase. In 
the vicinity where the three phases meet, the thermal 
signatures of the Photon to Higgs transition change dra- 
matically, see Fig. [IT] We interpret this as the transition 
becoming first order for JC < 0.2, which is also expected 
on general grounds due to increasing fluctuations of the 
molecular field, see Appendix [XJ This interpretation im- 
plies a tricritical point around JC ~ 0.2 separating the 
two regimes. 

It is difficult to characterize firmly the region of con- 
tinuous transitions in the model I. At each point, scaling 
apparently works, but the extracted exponents drift as 
we move along the phase boundary, with both v and r\ 
decreasing towards 0.5 and respectively as we approach 
JC = 0.2. We attribute this to crossovers near the tricriti- 



cal point. In this picture, the second-order region is bor- 
dered by the 0(4) transition at JC = 00 and the tricritical 
point near JC = 0.2. The RG flows near the correspond- 
ing fixed points are causing a wealth of crossovers. The 
absence of full leverage over such phenomena is the main 
source of our difficulties in this model. It is worth repeat- 
ing that in this picture the Higgs transition is controlled 
by a new fixed point, and there is no contradiction here 
since the RG flows are likely occurring in a more complex 
space than just the JC — J parameter space of the spe- 
cific model. By studying the dual stiffness and comparing 
with the NCCP model, we know that our systems are 
sufficiently away from the JC = 00 limit, and the appar- 
ently continuous transition with large v is most naturally 
explained as the second-order Higgs universality. 

To alleviate the difficulties associated with the 
crossovers near the putative tricritical point, we also con- 
sider a modification of the original NCCP 1 model by 
adding short-range antiferromagnetic interactions that 
hinder the formation of the Molecular phase. In fact, 
in the specific NCCP 1 model II, the Molecular phase is 
eliminated completely, allowing one to focus solely on the 
Photon to Higgs transition. We find that the transition is 
continuous along the whole phase boundary, see Fig. [TJ;. 
The scalings are consistent throughout, giving the large 
correlation length exponent. The transition line connects 
smoothly with the J = 00 limit, which can be viewed as 
a description in terms of the Abrikosov-Nielsen-Olesen 
vortices that have nontrivial internal structure. We are 
able to realize this structure and study the transition in 
Monte Carlo using local matter degrees of freedom, and 
this numerics gives our best estimates of the critical in- 
dices v and 77, Eqs. (|2~4"|) and (|2"5|) . For the future, it would 
be useful if one could realize such a nontrivial loop system 
in a manner that would eliminate critical slowing down, 
as is possible for the ANO vortices of the one-component 
model, which are short-range interacting loops with no 
internal structure. 

Let us conclude with more questions that one would 
like to study. First of all, we mention possible ap- 
plications to unusual quantum critical phenomena in 
magnetsj^^ Our correlation length and magnetization 
exponents are broadly in agreement with those found 
in the study of the continuous Neel to Valence Bond 
Solid transition in an SU(2) symmetric spin-1/2 system 
in Refs.0@. (It would also be useful to compare universal 
amplitudes as discussed in Ref. [H.) Motivated by this 
application, it would be useful to find the scaling dimen- 
sion of a monopole insertion operator in the gauge theory 
Eq. (UJ , since this corresponds^ to the VBS order param- 
eter measured in Ref££. Our J — 00 formulation of the 
model II, where inserting a monopole (antimonopole) cor- 
responds to a source (sink) for the discrete -B-fluxes, may 
be well suited for this. It would also be useful to find the 
scaling dimension of multiple monopole insertions. For 
example, one would like to check whether quadrupled 
monopolcs are irrelevant as proposed in the deconfincd 
criticality scenario for the Neel - VBS transition on the 
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square lattice and as needed for the non-compact gauge 
theory Eq. (Q]) to be applicable to this transition^ 

It would be also useful to know the effect of various de- 
formations away from the SU(2)-invariant case. For ex- 
ample, recently Alet et alZ studied a Coulomb to Valence 
Bond Solid transition in a classical dimer model on a 
3D cubic lattice. This can be (approximately) mapped^ 
to Eq. JT]) plus cubic anisotropy terms, which, however, 
appear only with more derivatives/fields. Alet et alZ 
found a continuous transition but with different critical 
indices from ours. Several possible explanations are given 
in Ref. 0, and more studies would be worthwhile. 

Finally, we would also like to revisit the U(l)xU(l) 
variant using the insights learned here. 
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Molecular phase, m^v becomes smaller leading to a sign 
change for the |\E'| 4 term. By the familiar mechanism 
in the Landau theory, the transition becomes first order 
due to such short-distance energetics near the molecular 
phase. 
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FIG. 14: Phase diagram of the Landau mean field theory 
Eq. (|A1|) in the vicinity where the three phases meet. The 
bold lines mark where the transitions become first order in the 
vicinity of the fork. We chose the axes so that m$, mjv > in 
the lower left corner, which gives the locations of the phases 
roughly similar to that in Fig. [Tp. 



APPENDIX A: MEAN FIELD ANALYSIS NEAR 
THE MOLECULAR PHASE 

Here we summarize the mean field argument that in 
the two-component system, the Photon to Higgs transi- 
tion becomes first order in the vicinity of the Molecu- 
lar phased We introduce an order parameter N for the 
molecular phase and write a schematic Landau theory 
that contains all three phases: 

S[9,rf] = m v ,|*| 2 + ^|*| 4 + m w iV 2 + ^(A? 2 ) 2 
- gN ■ * t ct* . (Al) 

Setting = un — g — 1 and minimizing this action, 
we obtain the phase diagram in Fig. 1141 Away from the 
fork all transitions are continuous, while near the fork 
the Photon to Higgs and Molecular to Higgs transitions 
become first order. Let us see how this happens for the 
Photon to Higgs transition. For large positive mjv, we 
can "integrate out" the N field and obtain the renormal- 
ization uy, — > — g 2 /(2mjv)- When we approach the 



APPENDIX B: EVALUATION OF THE HELICITY 
MODULUS T IN THE NCCP 1 MODELS 



The helicity modulus is defined from the dependence 
of the free energy on the twist in the boundary con- 
ditions. As is common when calculating stiffnesses in 
Monte Carlo, the twist angle can be distributed into small 
twists throughout the system, and the helicity modulus 
can be then obtained by measuring appropriate deriva- 
tive expressions during simulations with periodic bound- 
ary conditions. In this approach, the expressions to be 
evaluated depend on the detailed interactions and are 
different in the NCCP 1 model I and model II cases. 

Specifically, we consider twisting the vector ft = 
(711,71,2,71,3) in the n\ — n<i plane, which is achieved by 
the substitution {z r ^z r {) — ► {z r ^e llr ^, z r je _I7r ") in the 
action for the twist imposed in the lattice direction jl. 
Carrying through the derivations, we can summarize the 
expressions to be measured in the two models as 
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where we introduced two functions f(x) and g{x). In the 
model I case, f(x) = g(x) = 0. In the model II case, they 
are 



Ox) 



h{x) ^ 
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where Jo(x) and are the modified Bessel functions; 

the corresponding terms in the expression for T arise 
from the additional interactions defining the model II, 
Eq. ©. 
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